IU-MSTP/76 
Dec, 2008 



Magnetic translation symmetry on the lattice 



Ken-ichi Sekiguchi 1 , Tomohiro Okamoto 1 and Takanori Fujiwara 2 



1 Graduate School of Science and Engineering, Ibaraki University, Mito 310-8512, Japan 
2 Department of Physics, Ibaraki University, Mito 310-8512, Japan 



Abstract 

Magnetic translation symmetry on a finite periodic square lattice is investigated for an arbi- 
trary uniform magnetic field in arbitrary dimensions. It can be used to classify eigenvectors 
of the Hamiltonian. The system can be converted to another system of half or lower dimen- 
sions. A higher dimensional generalization of Harper equation is obtained for tight-binding 
systems. 



1 Introduction 



Systems of charged particles interacting with uniform magnetic fields have served as testing 
grounds of various ideas in quantum field theory and condensed matter physics. Uniform 
magnetic fields are simple and in many cases exact solutions are available while keeping 
some essential aspects of the interactions with the gauge fields. A remarkable feature of 
these system is the appearance of magnetic translation symmetry [T]. 

In the presence of magnetic fields the gauge potential depends on the coordinates. The 
change in the gauge potential due to a translation can be compensated by a suitable gauge 
transformation. Arbitrary translations become a symmetry of the Hamiltonian in infinite 
euclidean spaces and each eigenstate is infinitely degenerated. In finite periodic spaces or 
on compact tori the boundary conditions for the wave functions become twisted. Because of 
this only a discrete finite subgroup of the whole translations survives as magnetic translation 
symmetry. The degeneracy of states is a finte constant depending on the magnetic field. 

On finite periodic lattices exact solutions are not available and translation symmetries are 
rather restricted. The spectra are far richer than those of the continuum [2]. Nevertheless, 
it is possible to define magnetic translations, which enable us to qualitatively understand 
the richness of the spectra. In this paper we investigate magnetic translation symmetry on 
the lattice in full detail for arbitrary uniform magnetic fields in arbitrary higher dimensions. 
The key idea is to introduce oblique lattice where the magnetic field takes a block-diagonal 
form composed of 2 x 2 anti-symmetric matrices. According to the rank of magnetic field, 
the generators of magnetic translations can be classified into pairs satisfying two dimensional 
magnetic translation algebra and unpaired commuting generators. This enables us to identify 
maximal commuting subset of the generators. Representation theoretic approach [3] can be 
used to reduce Hamiltonians into block- diagonal forms without relying on its concrete form. 
The reduced Hamiltonian describes a system of half or lower than half dimensions. This 
phenomenon has been observed for systems of electrons in periodic potentials [I]. Magnetic 
translation symmetry has been investigated in four dimensions in the context of chiral lattice 
fermions [5] . Reduction from two to one dimension and the spectra of hermitian Wilson-Dirac 
operator have also been discussed in Ref. [6]. 

This paper is organized as follows. In the next section we show that periodic boundary 
conditions on finite lattice can be converted to twisted periodicity on infinite lattice. In 
Sect. [3J we examine magnetic translations on the lattice. Constraints on the translation 
vectors are solved. To illustrate the approach we argue magnetic translation symmetry in 
two dimensions. This is done in Sect. HI We introduce oblique lattice coordinates in Sect. 
[5] and give a generalization to arbitrary dimensions. Sect. [6] is devoted to summary and 
discussion. 

2 Twisted boundary conditions 

Let us consider a system of a charged particle interacting with a uniform magnetic field of an 
abelian gauge theory onarf dimensional periodic lattice of a size L d , where L is a positive 
integer. We take the lattice constant a = 1. Any uniform magnetic field F^ u = —F Vjx 
(//, v — 1, • • • , d) on the lattice can be written as 



F, 



2ixm 




L 2 



2 



where m^ v are integers. These classify topological sectors of the lattice gauge fields [3 [8]. 
The link variables U^(x) giving rise to the uniform magnetic field can be chosen as [7] 



exp 



2m 



,L-1 



E2iri ^-^ 



(2.2) 



where < x^ < L (pi — 1, • • • , d) stand for periodic lattice coordinates satisfying x M + L = x^ 
and 6^ are real constants. These satisfy periodic boundary conditions 



U^x + LP) 



(2.3) 



where v is the unit vector along the z/-th lattice axis. The magnetic fields (12 .11) are related 
to the plaquette variables by 

JFuu 



U^(x)U u (x + fi)UUx + v)Ul{x) 



(2.4) 



We denote the Hamiltonian or the lattice Dirac operator of the system by H and consider 
the eigenvalue problem 



Hil>(x) 



\ip(x). 



(2.5) 



Concrete form of H is not necessary to argue magnetic translation symmetries. We only 
assume that H only depends on the lattice coordinates through the link variables. The wave 
function ip{x) is subject to the the periodic boundary conditions 



ip(x + Lfi) = ip(x). 



(2.6) 



In lattice gauge theory it is customary to adopt periodic boundary conditions. The lat- 
tice is supposed to be a regularization of the infinite continuum and the boundary conditions 
become irrelevant in the infinite volume limit. In the present case the lattice should be un- 
derstood as a regularization of a finite volume flat torus and the periodic boundary condition 
is considered to be legitimate. In the continuum, however, neither parallel transporters, the 
continuum analog of lattice link variables, nor wave functions can be periodic on tori in the 
presence of a net magnetic flux. 

This apparent mismatch between the continuum and the lattice can be resolved by noting 
that the periodic link variables (12.21) become singular in the classical continuum limit a — > 0. 
The singularities can be removed by a gauge transformation 



A (x) 



exp 



IL 



T 



(2.7) 



where [c] stands for the integer part of c, i.e., [c] = n if n < c < n + 1 for some integer n. 
The Aq(x) removes the troublesome parts of ( 12.21) and yields 

Ufa) = A (x)Ufa)Al(x + fi) 
2ni 



exp 



exp 



L 2 
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This link variable has a smooth classical continuum limit and (12. 4p becomes apparent in 
this gauge. The cost we must pay is the simple periodicity of the link variable. The gauge 
transformation (12.71) is not periodic under the shift x — > x + L(l but satisfies 



A (x + Lfl) = Al(x)A (x). 




It changes the periodic boundary conditions to twisted ones. The gauge transformed wave 
functions 



A (x)if)(x) 



(2.10) 



are subject to 



^(x + Lft) = A a (x)ip a (x). 



(2.11) 



The link variable (I2.8p corresponds to axial gauge in the continuum. We can also work 
in the symmetric gauge by carrying out a further gauge transformation by 



A sa (x) 



exp 



7) / . " puV^^V 



p.<v 



;2.i2) 



The link variable in the symmetric gauge is given by 



XPJx) = exp 



(2.13) 



In this gauge the link variables and the wave functions satisfy the boundary conditions 

U;(x + LO)=Al(x)U;(x)A s J(x + tt, r(x + Lfi) = A*(x)V> s (x), (2.14) 
where A s ^(x) is defined by 



exp 



%L 



Tf, 



2 ^ 



(2.15) 



This gauge is suitable for analyzing magnetic translation symmetry in the next section. 

The constant phases 6 M in the link variables do not affect the magnetic field F^. In the 
continuum we can remove them by a suitable choice of the coordinate origin if det F^ v ^ 
0. The spectrum of the Hamiltonian does not depend on them. This is not the case for 
det F^ u = as is observed in odd dimensions. In lattice gauge theory we can also remove 
bfj, from the link variables by a further twisting of the boundary conditions. The eigenvalues 
depend on even if det F^ u ^ is satisfied. In this section we have not carried out this 
twisting, leaving the constant term intact. 
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3 Magnetic translation symmetry on the lattice 



In the continuum the Hamiltonian of a charged particle in a uniform magnetic field is not 
invariant under arbitrary translations. This is due to the coordinate dependence of the gauge 
potential. A discrete subgroup of the translations known as the magnetic translation p], 
however, survives. It is a suitable combination of translations and gauge transformations. 
We expect that a similar situation also occurs on the lattice. In this section we pursue the 
conditions for the magnetic translations in lattice gauge theory. 

We work in the symmetric gauge (12.131) and consider a shift of the lattice coordinates 
x — > x + £, where £ is an arbitrary integer vector. It is easy to verify that the link variables 
U^(x) and U*(x + £) are related by 



where fi s e (x) is given by 



(3.1) 



n s Ax) 



exp 



2 Fv.v^n x v 



(3.2) 



The relation (13.1 ft is analogous to a gauge transformation. In general we cannot regard 
U^(x + 1) as a link variable associated with the link (x,fi) since it does not satisfy the 
boundary conditions (I2.14p . However, eq. (13.11) suggests a transformation TJ» : ip s — ► ijf t 
defined by 



i)\{x) 



ft? I 



x)4> s {x + £). 



(3.3) 



If i/j s (x) is an eigenvector of the Hamiltonian and ip\{x) satisfies the twisted boundary condi- 
tions ( 12.141) . then ip}(x) is also an eigenvector belonging to the same eigenvalue with ip s (x). 
In other words the transformation ( 13.31) is a symmetry of the Hamiltonian. We call this as 
magnetic translation symmetry on the lattice. 

The requirement that (13. 3p satisfies the boundary conditions (12.141) leads to the consis- 
tency conditions for f2|(x) 



A s (x + £)n s e (x)=n s e (x + mA s (x). 



(3.4) 



This can be seen by noting that a translation x — > x + L£l + £ can be achieved in two different 
ways x — > x + Lfi — >• x + Lfi + £ and x— > x + £ — > x + £ + Lfi. The conditions (13.41) can be 



stated in terms of as 



2_, m ^u£u = mod L. 



(3.5) 



One can also arrive at the same result by working in other gauges. In the axial gauge f 1 2 . 8 1) 
magnetic translations are implemented by 



exp 



^ ^ ^ P\xv£ [1%P 



(3.6) 
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We will use this later. 

We now turn to analyzing the constraints (I3.5p . We can block-diagonalize into 2x2 
anti-symmetric integer matrices as in the continuum [HJ El ED] 



m 



peri 



(3.7) 



p.a 



where C = {C^) is an integer matrix with det£ = 1 and v = [y^v) takes the form 

/ vi \ 
-v x 

v m 
-v m 



V 



(3.8) 



0/ 



The set of integers v p (p = 1, • • • , m) with 2m being the rank of (m Mi ,) can be chosen so that 
u p divides v p+ \ for p — 1, • • • , m — 1. 

We write £ _1 in terms of d integer vectors M p , N p and K\ (a = 1, • • • , m, I = 2m + 
1 , • • ■ , d) as 



CT l = (Mi, Ni, ■ ■ ■ , M m , N m , K 2m+1 , • ■ ■ , K d ) T 
We also introduce the dual integer vectors M* p , N* p and K* 1 by 



£ = (M*\ N*\ ■ ■ ■ , M* m , N* m , K* 2m+1 , ■ ■ ■ ,K* 



Then f !3.7p can be written as 



V ) 



(3.9) 



(3.10) 



(3.11) 



The set of d integer vectors M, N and K also generates the original lattice. We can expand 
an arbitrary integer vector £ as 



d— 2m 



p=l 



(3.12) 



i=i 



where £ P M , £ P N and i l K are all integers. The constraints (13. 5p for the magnetic translations 
can be simplified as 



v p £ p M = v p £ p N = mod L, 



(3.13) 



whereas £ l K can be arbitrary. 

To solve the conditions (13.131) we introduce a set of integers r p , s p > and n p by 



'pi 1 



(3.14) 
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where n p and s p are mutually prime. Then £ P M = js p and £ P N = ks p satisfy (13 . 13[) for arbitrary 
integers j, k. We thus find three basic vectors generating magnetic translations 

£ ( p M) = s p M p , £^ = s p N p , &* ) = K l . (3.15) 

In the following sections we need another set of integers r' p , s' p > and n' p defined by 

n'% = n p , r' p s' p = r p , (3.16) 

where n' p and r' p are mutually prime. 



4 Magnetic translations on 2d lattice 

To illustrate magnetic translations on the lattice we consider a two dimensional system with 
a magnetic field = 2iti>/L 2 . This is already of a standard form (13. 7p and £ is taken to 
be the identity matrix. The integer vectors Mi, iV\ and their dual are given by 

Mi = M* 1 = , Ni = N* 1 = 

As in (I3.14p . we have a unique set of integers n = n 1 , r = r\ and s = s 1 for a given 
L and v — v\. We also define another set of integers n' = n n , r' = r[ and s' = s' 1 by 
(I3.16p . We thus find two independent basic magnetic translations T x = T sMl and T y = T sNl 
corresponding to x — > x + sM\ and x — > x + sNi, respectively. In the axial gauge these act 
on the wave functions by 

T x ip(x,y) = e 3 ^ip(x + s,y), T y tp(x, y) = tp(x, y + s), (4.1) 

where use has been made of (13.61) . We have suppressed the label "a" standing for the axial 
gauge. As can be verified from (14.11) . these operator satisfy 

T x T y = e~^T y T x . (4.2) 
This gives T x T y = T y T x ' . We can find simultaneous eigenvectors of if, T x and T y . The 

/ 2iirq 

eigenvalues of T£ can be written as e~^~ by an integer q (0 < q < s ) since T x satisfies 
(T; Y = T x r = 1. Similarly, we can write the eigenvalue of T y as e r with an integer 
< p < r. Let us denote a simultaneous eigenvectors by i[> p ,q, then we have 

Hip P , q (x,y) 

Ty1pp,q(X } y) 

The magnetic translation T™ (m = 0, 1, • ■ ■ , r' — 1) maps an eigenstate of T y belonging to an 

eigenvalue e r to other eigenstate with an eigenvalue e r .In other words, it changes 
the label p mod r to p + ran mod r. Since p + r'n = p mod r, there are r' degenerate 
eigenvectors belonging to the eigenvalue X P) g. We thus find that the eigenvalue X p>q depends 
only on p mod s' with respect to p. 




\p, q ipp, q {x,y), 
^ P ,q(x,y), 



2 1. it q 



2iTvp 



e r ip p Jx,y). 



(4.3) 
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The second and the third of (14. 3 j) together with (14. ip imply that the wave functions at 
(x, y) and (x + sr', y) or (x, y + s) are related by 

. 'Hivn 1 y ,2iivq 2iirp 

tpp, q {x + sr ,y) = e ~ + ^ip p ^{x, y), ip P , q [x,y + s) = e~tfj Pjq (x,y). (4.4) 
Noting the twisted periodicity in y, we can expand ip p , q (x, y) in Fourier series as 



1 s_1 

^ p , q {x,y) = -y=J2vp^j( x ) e 



2i7r(p+jn'r)y 
L 



(4.5) 



j=0 

where ip p , q ;j(x) must satisfy 

¥ P , q]j +s(x) = (p Pl q;j{x), ¥>p, q -A x + sr ') = e ~^Vp,«i+i(aO- (4.6) 

We see that the original L 2 components of ip P ,q{x,y) (0 < x, y < L) can be expressed in 
terms s 2 r' values of a one dimensional wave function ip m (x) = <p p , q ;o{x) (0 < x < s 2 r'). It 
satisfies the twisted periodicity 

(p P , q (x + s 2 r') = e^cpp !q (x) (4.7) 

as one can see from (14.61) . The magnetic translations put no further restrictions on (p pq (x). 

Our arguments so far does not depend on the detailed form of the Hamiltonian H. We 
can block-diagonalize it into r' matrices whatever form it is. This leads us to the Hamiltonian 
H p -j defined by 

_L ^ — T 2i7v(p J rjn r)y 

H^p, q {x,y) = —/=2^ e Z n p;j¥ P ,q-d( x )- ( 4 -8) 

It acts on the one dimensional system of a size s 2 r' satisfying the twisted boundary condition 
(14.71) . Since the eigenvalues of H p -j do not depend on j, it is enough to consider the case 
j = as mentioned above. We thus arrive at the eigenvalue equation for (p PA (x) 

np<Pp,q( X ) = VsPmO*)' (4-9) 

where TC P stands for TC p -o- 

For tight-binding Hamiltonian ip Piq (x,y) satisfies 



e ibx *p P ,q(x + l,y) + e~ ih *%, q {x -l,y) 

2%Ttn x i/L , _ Znrn x ■ l 

+e ^ +4b ^ P , g (x, V + 1) + e W lb ^ m {x, y - 1) = \p, q ip p , q (x, y). (4.10) 



The equations for ip Pjq (x) can be found by inserting (14.51) into the expression. The constant 
phase factors e ±lbx can be eliminated if we introduce 

<p m {*) = <P"<pM- ( 4 - n ) 

We thus obtain 

f / 2irn'x 2itr> \ 1 

<p P)9 {x + 1) + y p , q (x - 1) + I 2 cos f —^j- + + b y 1 - X p>q J. £ M (a;) = 0. (4.12) 
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-3 - 




7t/6 7t/3 Jt/2 

Parameter b r 



Figure 1: Eigenvalues of tight-binding Hamiltonian in two dimensions are plotted for < 
b x < 7r/2 and b y = 0. We take L — 6, v — 9. 



This is referred to as the Harper equation [TT] in condensed matter physics. The periodicity 
(14.71) under the shift x — > x + s 2 r' is further twisted by (14. lip and is given by 

<p M (x + s 2 r') = e^ +ib ^ r '<f Piq (x). (4.13) 

The plot of the spectrum as a function of a = n'/s 2 r', i.e., magnetic flux per plaquette, 
is known as the butterfly diagram [2]. Unlike the continuum theories eigenvalues depends 
on the parameters b x and b y . The twisted periodicity (14. 13j) is invariant under the shift 
b x ~~ > b x + 2n / s 2 r' . This implies that the eigenvalues are periodic in b x and also in b y by 
rotation symmetry with a period 2n/s 2 r'. See Fig{TJ 



5 Extension to higher dimensions 

We have shown that the magnetic translation symmetry can be used to constrain the wave 
functions on a two dimensional square lattice. We now extend this to higher dimensions. 

Magnetic translation symmetry is apparent for special types of magnetic fields like (13. 8p 
and the analysis for the two dimensional system in the previous section is applicable imme- 
diately. This is not the case for a general uniform magnetic field. Fortunately, it is always 
possible to transform F^ v into the block-diagonal form (13. 8p as mentioned in Sect. [31 In 
the continuum eigenvalue problems on magnetized tori can be solved by using this approach 
|12j . In this section we use the notation of Ref. [TO] . 

We introduce oblique lattice coordinates £ p , rf (p = 1, ■ ■ • , m, / = 2m + 1, • • • , d) by 

x = j2( M P e + nx) + E K u l > t 5 - 1 ) 

v i 
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where M p , N p and Ki are oblique lattice vectors defined by (13. 9p . The oblique lattice 
coordinates are all integers and have one-to-one correspondence with x as can be seen from 



X 



E A ", 



* l x 



(5.2) 



In particular the unit translations x —>■ x ± jl correspond to the following shifts 



f p -> £ p ± M* p , 



if ± 



*/ 



(5.3) 



Conversely, the translation £, 77, x — > £ ± e p , 77, x can be realized by the shift x — > x ± M p , 
where e p is the p-th unit vector with (£ + e p ) 9 = £ g + <5 p>g . Similar things hold true for the 
unit translations in rj or in 

The link variables in the symmetric gauge (I2.13P can be written as 



Kip) 



exp 



cxp 



—I 



(5.4) 



where n' p , r p and s p are given by (13.141) and (I3.16p . To argue magnetic translations it is 
more convenient to work in axial gauge on the oblique lattice. This can be achieved by the 
following gauge transformation 



U"{x) = A as {x)U s Jx)A as \x + p) = exp 



y ^£Ln? (e+ -m*A + ib, 



with A as (x) = exp 



—1 



E 



7T77, 



(s p ) 2 r p 



7 ev 



(5.5) 



In this gauge the twisted boundary conditions for the wave functions are given by 



with A°(rr) = exp 



E 



2vm p 
sp 



K p v p + -zN? 



(5.6) 



We have suppressed the label a indicating the axial gauge for the wave function. 

We can also find the periodicity of the wave function under shifts of the oblique lattice 
coordinates by L by a repeated use of (15. 6p as 



+ Le p ,r},x) 



AN) 



■^(t,V + Le p ,x) 
^(£,ViX + Lei) 

where e's are integers defined by 



(5.7) 



jl<V 



jl<V 
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We see that the wave function is either periodic or anti-periodic in rf and x ■ 

We now turn to magnetic translation of the wave function. In the axial gauge (15.51) the 
link variables on the sites x and x + i are related by 

m(x + i) = n«(x)u a Jx)nf(x + fi) 



with VI i {x) = exp 



V 



where £ is given by (13.121) . This leads to magnetic translation of the wave function 

T*l>{t V, X) = Vif(x)^ + £ M , V + £n,X + tie). (5.9) 

Let us denote the generators of magnetic translations corresponding to the three basic trans- 
lations (13.151) by Tp M \ Tp N ^ and T^ K \ respectively. Then (15. 9p yields 

2mn' p 



s p r 



— T 



T p (M V(e,^,x) = exp 

T i V, x) = ^(6, V, X + ez), 

The generators of magnetic translations satisfy 

2inn'P 



(5.10) 



(5.11) 
(5.12) 



-p ~p - p p 

and all other combinations are commutative. They also satisfy 

^ T {M)y' p y' p = ? (TWfv = e ™4 N) ? (T} K) ) L = , 

as can be seen from ( 15.7ft and ( 15.101) . These are a higher dimensional generalization of ( 14. ip 
and (jOl . 

We may choose a set of commuting operators H, (Tp ) r 'p, T p N \ 7] as in two dimensions 
and consider eigenvectors of these operators 

Hlpk,p,q(€,V,X) = h,p,q^k,p,q(€,V,X), 

'2m ( 1 



(TW)^ fcjM (e,r7,x)=exp 
T p (Ar V^, g (e^,x)=exp 
Ti^k, P ,q(C,V,X) = exp 



2ztt 

2«7T 



(M) 



(K) 



*Pk,p,q(^V,X), 
tpk,p,q(^V,X), 



(5.13) 



where < q p < s' , < p p < r p and < k\ < L are integers. These together with (15. 101) 
implies 



^k, P ,q{i + s p r'e p , 7], x) = exp 



2iim' p „ 2ztt 



S (A0 



ipk, P ,q(£, V + «%, x) = exp 
V>fc, P ,g(£, »7, X + ej) = exp 



— (p + - e w 



2m 



^k,p, q {i, , n,x)i 

*Pk,p,q(LV,X)- 



(5.14) 
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The last two of these expressions enable us to write 

A, r „((,V,x) - £» W ({)e*fe(^^W" , )«'^(''*«' n M > (5.15) 

Vs 1 • • • s m ~ 

3 

where the sum is taken over integers < j p < s p (p = 1, • • • , m). The Fourier coefficients 
(Pkj>,q;j{Q m ust satisfy the following periodicity 

<Pk,p,q;j+sre p (€) = <Pk,p,q;j(€), ¥k,p,q;M + ^Cp) = C S " V " ' " > fk,p,q;j+e p {0 (5.16) 

The first comes from the definition of <Pk,p,q;j{0 arid the second from the first of (I5.14p . We 
see that the original L d components of the wave function ipk, P , q (x) (0 < x M < L, \i = 1, ■ ■ • ,d) 
can be obtained from the (s 1 ) 2 r^ • • • (s m ) 2 r' m components of <fk,p,q(£) = V ? fe,p,g;o(0 (0 — £ P < 
{s p ) 2 r' p , p = 1, • • • ,m). It must satisfy the twisted boundary conditions 

Vk,p,q{t + ^ P ?r' p ep) = e^T^ P ^WO (5.17) 

and can be considered as a wave functions of some m dimensional system. The Hamiltonian 
7~tk, P ;j of the dimensionally reduced system is found by inserting (15.151) into the first of (15.131) . 
i.e., 

V s • • • s m 

1 

(5.18) 

It is enough to consider the case j = 0. We thus obtain 

7~Lk,p l Pk,p,q{.0 = ^k,p,q^k,p,q(0^ (5.19) 

where we have defined TCk, p = 7~Ck,pfi- This together with (15.171) determines the wave func- 
tions <fk,p,q, or equivalently <Pk,p,q;j-> an d the eigenvalues \k, P ,q- They are defined on the m 
dimensional oblique lattice of a size (s 1 ) 2 ^ x • • • x (s m ) 2 r' m . We can then reconstruct ipk, P , q 
by (15.151) and arrive at r\ • • -r^s^ • • • s' m L d ~ 2m wave functions forming a complete set. 

For any fixed k and q there are r± ■ ■ -r m wave functions i>k,p,q- They can be decomposed 
into s' x • • • s' m sets of wave functions degenerated by the magnetic translation symmetry. For 
a given set of labels p q mod s' (q = 1, • • • , m) there are exactly r[ - ■ ■ r' m degenerate wave 
functions given by 

{T (M ))hl . . . ^lO)^. (0<h q <r' q ,q=l,...,m) (5.20) 
We now apply the formalism developed so far to the tight-binding system described by 

d 

^2{ U ^( X )^k,p,q(x + P) + U^(X - fl)^ k ,p,q(x ~ p)} = \ k ,p,q1pk,p,q{%), (5.21) 
ft=l 

where the link variables are given by (15. 5ft . This can be converted to the form (15.191) . To 
have a compact expression we introduce /3's by the expansion 

K = Y,^ I)M 7 + ft" )N ?) + E A W <> (5-22) 

q I 
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and define p, q and k by 
Then (15.211) is reduced to 

+6 -¥E,^{?H«?W-¥s^^ a _ M ;)} = 

(5.24) 

The phase factor containing can be removed by introducing 0k,p,q as 

WO = e'^'^^VWfl- (5-25) 

We finally obtain 



+e -^E 8 ^H^(^-|M-) + p 9 }-^ E! ^^^ (e _ M *) J 



^k,p,q^Pk,p,q(jC) j 



with ^ fejM (£ + (s p ) r p e p ) = e s * P <Pk, P , q (0 (5.26) 

This is a higher dimensional extension of the Harper equation (14.121) . The twisted pe- 
riodicity can be derived from (15.171) and (I5.25p . It is periodic under the shift — > 
p(M) _|_ 2wep/(s p ) 2 r' p . This implies that the spectrum is symmetric under the same shift as 
has been observed in two dimensions. 



6 Summary and discussion 

We have investigated magnetic translation symmetries on finite periodic lattices in arbitrary 
dimensions. We have shown that any lattice system with uniform background magnetic 
field possesses the symmetry. In continuum theories magnetic flux tensor determines 
a unique magnetic translation group. In lattice theories, however, it depends on how the 
lattice size L is taken. This is due to the simple fact that only translations by a divisor of L 
is allowed on the lattice whereas u p is not necessary a divisor of L. The symmetry changes 
rapidly as one varies the magnetic field. Then the spectrum of H is very sensitive with 
respect to the applied magnetic field. These qualitative differences between the continuum 
and the lattice disappear in the continuum limit. We can always take the limit keeping any 
prescribed magnetic translation symmetry. 

In the extreme case that there is no magnetic translation shorter than the period of 
the lattice, eigenvectors cannot be constrained by the magnetic translation symmetries. We 
cannot reduce the number of unknown variables in solving the eigenvalue problem. However, 
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it is always possible to dimensionally reduce the system. An immediate consequence of this 
is that any tight-binding system can be formulated by the generalized Harper equations. 
Dimensionally reduced systems may be more tractable than the original equations and be 
beneficial to further studies. 

As an application of magnetic translation symmetry we can compute index of overlap 
Dirac operator for abelian gauge background. It is related to spectral asymmetry of hermitian 
Wilson-Dirac operator. Due to the topological invariance of the index it suffices to compute 
the spectral asymmetry for a uniform magnetic field belonging to the same topological sec- 
tor. Since the degeneracy of each eigenvector is easily seen from the magnetic translation 
symmetry, we have only to compute the spectral asymmetry of the reduced Hamiltonian. 
This has been carried out in two dimensions. Extension to higher dimensions is certainly 
interesting. We will argue the magnetic translation symmetry in lattice fermion systems 
elsewhere. 
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